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The pulse profile and spin evolution of the accreting pulsar in 
Terzan 5, IGR J17480-2446, during its 2010 outburst 
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ABSTRACT 

We analyse the spectral and pulse properties of the 11 Hz transient accreting pulsar, IGR 
J17480-2446, in the globular cluster Terzan 5, considering all the available Rossi X-ray Timing 
Explorer, Swift and INTEGRAL observations performed during the outburst shown between 
October and November, 2010. 

By measuring the pulse phase evolution we conclude that the NS spun up at an average 
rate of < ;> >= 1.48(2) x 10^^^ Hz s^^, compatible with the accretion of the Keplerian 
angular mom entum of ma t ter at the inner disc boundary. This confirms the trend previously 
observed by Papitto et al.l (1201 lb . who considered only the first few weeks of the outburst. 
Similar to other accreting pulsars, the stability of the pulse phases determined by using the 
second harmonic component is higher than that of the phases based on the fundamental fre- 
quency. Under the assumption that the second harmonic is a good tracer of the neutron star 
spin frequency, we successfully model its evolution in terms of a luminosity dependent ac- 
cretion torque. If the NS accretes the specific Keplerian angular momentum of the in-flowing 
matter, we estimate the inner disc radius to lie between 47 and 93 km when the luminosity 
attains its peak value. Smaller values are obtained if the interaction between the magnetic field 
lines and the plasma in the disc is considered. 

The phase-averaged spectrum is described by thermal Comptonization of photons with 
energy of w 1 keV. A hard to soft state transition is observed during the outburst rise. The 
Comptonized spectrum evolves from a Comptonizing cloud at an electron temperature of « 
20 keV towards an optically denser cloud at kTf, « 3 keV. At the sam e time, the puls e am- 
plitude decreases from 27% to few per cent, as already noted by Papitto et al.l (1201 lb . and 
becomes strongly energy dependent. We discuss various possibilities to explain such a be- 
haviour, proposing that at large accretion luminosities a significant fraction of the in-falling 
matter is not channelled towards the magnetic poles, but rather accretes more evenly onto the 
NS surface. 
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1 INTRODUCTION 

The study of the frequency variations of the coherent signal emitted 
by an accreting pulsar is one of the foremost techniques to study the 
dynamical interaction between the rotating neutron star (NS in the 
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following) and the in-falling matter. This is particularly the case of 
disc-fed pulsars since the torque acting on the NS is expected to 
reflect the specific angular momentum of the disc plasma before it 
is captured by the pulsar magnetosphere, as well as the interaction 
between the disc pl asma and tho se magnetic field lines that close 
through it (see, e.g.. lGhoshll2007l and references therein). 

The process of disc angular momentum accretion is also es- 
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tablished as the driver to accelerate a NS until it reaches a spin 
period of few ms. According to the r ecycling scenario (see, e.g., 
iBhattacharva & van den Heuvellll99lh rotation-powered millisec- 
ond pulsars are spun up by a long-lasting phase of accretion of mass 
and of the relative angular momentum by a NS in a low-mass X-ray 
binary (LMXB). The discovery of coherent ms pulsations in the X- 
ray emission of a transient LMXB ( Wiinands & van der Klis 1998) 
provided a fundamental confiiTnation of this scenario, also open- 
ing the possibility of studying the frequency evolution of quickly 
rotating NS during the accretion phase. However, so far measur- 
ing the NS spin evolution while accreting of the known 14 ac- 
creting millisecond pulsars (AMSP) has been made often difficult 
by the intrinsic small torques exerted on the NS of these systems. 
This is partly due to the relatively short timescales of the accretion 
episodes, with outbursts lasting generally less than few weeks, and 
to the relatively small peak outbursts luminosities, which are not 
observed to exceed a level of few xlO"^'' erg s~^. To complicate 
further the task, timing noise was observed to affect the phase of the 
pulse profile (see, e.g.. lBurderi et al.l2006tlHartman et al.l2009r) . A 
spin-up at a rate of the order of that expected according to standard 
accretion theories could then be confirmed only in a couple o f cases 
JFalanga et alj2005l : lBurderi et al.l2007l;|PaDitto et al.l2008h . 

Here we present the case of the 1 1 Hz accreting pulsar, IGR 
J17480-2446, belonging to the globular cluster Terzan 5. The tem- 
poral evolution of the 1 1 Hz spin frequency during a subset of all 
the observations performed by th e Rossi X-Ray Timin g Explorer 
(RXTE) and Swift was studied bv lPapitto etaiTi201ll. Pll in the 
following), who found an average spin up rate compatible with ac- 
cretion of the Keplerian angular momentum of disc matter. The 
analysis of the Doppler shifts of the signal frequency allowed to 
measure the 2 1 .27 hr orbital period and to constrain the mass of the 
companion between 0.4 and 1.5 Mq (see Pll, and th e references 
there in). A candidate IR counterpart was proposed by iTesta et al.l 
( 1201 1) thanks to an observation performed by ESO-VLT/NAOS- 
CONICA while the source was in outbur st. This candidate was 
then identified bv lPatruno & Milond ( 1201 2h in a HST archival ob- 
servation performed when the source was in quiescence. From the 
presence of pulsations at the minimum and maximum luminosity 
Pll estimated the surface magnetic field of the NS to lie between 
a few times 10* and a few times 10^" G. These properties make 
this recently detected mildly-recycled pulsar, the only discovered 
so far in a LMXB to spin at a frequency between 10 and 100 Hz, a 
unique case study to directly observe the effects of disc torque on 
the NS through the measure of its spin evolution while accreting. 
In this paper, we analyse in detail the properties of the energy spec- 
trum (Sec. 13. il l and of the pulse profile (Sec. l3.2t of J17480 using 
all the available X-ray observations performed by RXTE, Swift, and 
INTEGRAL, focusing on the effects of the torque exerted on the NS 
by the disc accreting matter (Sec. l3.3l and l3.4) . The results obtained 
are then discussed in Sec.|4] 



2 OBSERVATIONS 

In this paper we consider all the observations of IGR 117480-2446 
(J 1 7480 in the following) performed by RXTE, Swift and INTE- 
GRAL during the intensive observational campaign w hich followed 
the d iscovery of the source on 2010 October 10.365 (IBordas et al.l 
l20ld : MJD 55479.365, all the epochs are given in UTC). The mon- 
itoring of the source ended at MID 55519.130, due to a solar oc- 
cultation of the sky region where this transient lies. Such obstruc- 
tion ceased only ~ 2 months later (MID 55584.562); since further 



RXTE observations failed to detect an excess above the background 
it has to be concluded that the 117480 outburst ended during the pe- 
riod of non- visibility. 

During the outburst the source showed a large number of X- 
ray bursts, whose spectral and temporal properties are compat- 
ible with thermonucl ear burning of a mix ed H/He layer on the 
NS surface (Motta et al.l20IlMChakrabortv & Bhattacharvval20Ill : 



iLinares et al.ll20I ih . Since we are here interested in the analysis of 
the source properties during the non-bursting emission, we discard 
~I0 s prior the onset of each of the bursts that could be signifi- 
cantly detected above the continuum emission, and a variable inter- 
val of length between 50 and 160 s after, depending on the burst de- 
cay timescale. During the analysis of observations taken in modes 
with low temporal resolution (e.g. Standard 2 modes of PCA aboard 
RXTE) the length of the time intervals excluded from the analysis 
was increased in order to equal the time-binning. 

2.1 RXTE 

In this paper we consid er data obtained by the Proportional Counter 
Array (PCA; 2-60 keV llahoda et al.l200(5h during the interval MID 
55482.008-55519.130 (Obsld. 95437). Spectra and response ma- 
trices were extracted by using standard tools of the Heasoft soft- 
ware package, version 6.11. The 'bright' model was used to sub- 
tract background emission, and only photons recorded by the top 
layer of the Proportional Counter Unit (PCU) 2 and falling in the 
2.5 - 30 keV energy range were retained in spectral fitting. Further, 
a systematic error of 0.5% was added to each spectral channejj. 
HEXTE data were not included in the analysis since the background 
dominates over source photons, and a number of 'line-like' resid- 
uals of calibration origin make the characterisation of the spectral 
continuum at high energies extremely problematic. In the temporal 
analysis we retained photons collected by all the PCUs of PCA in 
the 2-60 keV energy band, and encoded in fast configurations such 
as Good Xenon and Event, with temporal resolution of I and 122 
1J.S, respectively. 

2.2 Swift 

The Swift observations of 117480 covered the interval MID 
55479.198-55501.559 (Obs.Ids 31838, 31841, 437313, 437466) 
and in this paper we consider only data observed by the X-ray 
Telescope (XRT; 0.2-10 keV; i Gehrels et al."2004'). All the data 
were processed by using standard procedures (BuiTows et al. 2003) 
and the latest calibration files available at the time of the analysis 
(2011 September; caldb v. 20110725). Spectral analysis was per- 
formed on data collected both in window timing (WT) and photon 
counting (PC) mode (processed with the XRTPIPELINE v.0.12.6), 
while only WT data sampled every 1.7 ms, and falling in the 0.2- 
10 keV, were considered for the temporal analysis. Filtering and 
screening criteria were applied by using FTOOLS (Heasoft v.6.1 1). 
We extracted source and background light curves and spectra by 
selecting photons in the 1-10 keV energy band and event grades 
of and 0-12, respectively for the WT and PC mode. Grades 1 
and 2 in the WT data were filtered out in order to limit the effect 
of the recently reported redistribution problem on the source and 
background product (see the web-page http://www.swift.ac.uk for 
details). The only observation considered that was performed in PC 
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mode was affected by a strong pile- up, and correc ted according to 
tlie teclinique developed by Vaugha n et alj ( |2006[) . Given the rel- 
atively high count-rate of the source during the WT observations, 
we also checked if these data were affected by pile-up by extract- 
ing for each of them a source spectrum from a circular region in 
which we progressively excluded a larger and larger fraction of the 
inner central pixelo Fits to these spectra did not reveal any sig- 
nificant change in the model parameters, and thus we assumed that 
XRTAVT data were not affected by pile-up. Spectral channels were 
grouped to contain at least 50 counts each, as it was done for all 
the spectra produced by the various instruments considered in this 
paper, and a systematic error of 3% was added to each spectral bin. 



2.3 INTEGRAL 

INTEGRAL observations are commonly divided into "science 
windows" (SCWs), i.e., pointings with typical durations of ~2- 
3 ks. We considered all available SCWs that were performed in 
the direction of the source during the outburst. These comprised 
publicly available observations of the Galactic bulge (GB) region 
in satellite revolutions from 975 to 980 (Obsld. 0720001; MJD 
55479.365-55491.940), and the data collected during the ToO ob- 
servation of PKS 1830-211 in rev. 981 (Obsld. 0770005; MJD 
55494.778-55497.447). During the GB observations, J17480 was 
relat ively close to the centre of the IBIS/ISG RI telescope (20-150 
keV; lLebrun et al.l2003l ; IUbertini et al.l2003[) and thus it was simul- 
taneously observed also with the JEM-X telescope (4.5-25 keV; 
[Lund et al. 2003). During the observations performed in rev. 981, 
the source was always outside the JEM-X field of view. All INTE- 
GRAL data were analyse d using version 9.0 of t he OSA software 
distributed by the ISDC iCourvoisier et al.ll2003h . JEM-X and IS- 
GRI spectra were rebbinned in order to have 16, 32 or 64 energy 
channels in the different revolutions, depending on the available 
statistics of the data. 



3 ANALYSIS 

3.1 Spectral Analysis. 

The spectra of the 46 observations performed by PCA were 
modelled using a thermal Comptoniza tion model (nthcomp, 



mogelled using a tnermal L-omptoniza tion model (nthcomp, 
IZdziarskiet"aiTi 19961 ; IZvcki et all 1 19991) . with the addition of a 
Gaussian emission line centred at energies compatible with Iron K- 
a transition (6.4—6.97 ke V). Interstellar a bsorption was modelled 
using the TBAbs model jWilms et al.ll200 0), fixing the absorption 
colum n to the value determ ined thanks to a 30 ks Chandra observa- 
tion by JMneLetZJIloill), A^jj = 1.17x10^^ cm"^. The assumed 
model describes satisfactorily the observed spectra (Xr — 1; see the 
rightmost column of Tab. |3] where the parameters of the spectral 
continuum, fluxes in the 2.5-25 keV energy band, and luminosity 
extrapolated to the 0.05-150 keV band, are also reported for each 
spectrum). 

The nthcomp model describes a Comptonized spectrum in 
terms of the cloud electron temperature , fcTe, of the temperature 
of photons injected in the cloud, kTs and of an asymptotic power- 
law index, F, related to the cloud optical depth, r, by the relation 



r = 



9 1 

4^ T{l + T/3){kTe/m,c^) 



-,1/2 



1 
2' 



(1) 





10 




9 




8 


. — . 




'cfl 


7 




6 


o 


5 


to 


4 
3 




2 




1 







> 

g 


1.1 
0.9 



0.7 



> 25 



10 
5 




-I 1 r 




^>T^v^At^|W-^''^v.^tth^ 




5 10 15 20 25 30 35 40 

Time (days) since Tq (IVIJD 55479.81229) 

Figure 1. (a) X-ray luminosity of J17480 estimated by extrapolating 
the best-fitting spectral model, TBAbs* (nthcomp+gaussian) , to the 
0.05-150 keV energy band. Black points refer to the modelling of the ob- 
servations performed by PCA on-board RXTE, alone, red triangles to si- 
multaneous XRT-S'n'i/i' and PCA-RXTE observations, while green squares 
to simultaneous XRT-Swift, JEMX-INTEGRAL and ISGRl-INTEGRAL 
observations. The temperature of the soft input photons for Comptoniza- 
tion, the electron temperature of the Comptonizing region, and its optical 
depth, are plotted in the panels (b), (c) and (d), respectively. Plotted errors 
are given at the 90% confidence level. The solid line in panel (a) is a cubic 
spline approximation. 
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(see, e.g.. lLightman & Zdziarskilll987h . The optical depth was then 
estimated by inverting such a relation, and the relative uncertainty 
estimated by propagating the statistical errors on fcTe and F. The 
values of the parameters of the Comptonized spectrum measured 
across the outburst , kTs, kTe, and r are shown as black points 
in panels (b), (c), and (d) of Fig. [T] respectively. The 3-30 keV 
spectrum of the source is relatively soft if compared with AMSP 
(see, e.g.. IPoutaneal2006l , and references therein). With the excep- 
tion of the first RXTE observation (labelled as [B] in the top panel 
of Fig. [1] and spanning the interval MJD 55482.008-5482.043), 
the electron temperature of the Comptonizing medium is in fact 
relatively low ( ^ 3.5 keV) and the optical depth relatively large 
(r « 8-10). By imposing energy conservation between the cold 
(kTs ~ 0.9 keV) and the hot phase (kTe ~ 3 keV during most 
of the observations) producing the Comptonization spectrum, an 
estimate of the radius of the region that provides the soft pho- 
tons that are subsequently Comptonized in the hot cloud is oh- 
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Figure 2. Unfolded spectra accumulated during Swiff and INTE- 
GRAL observations covering the intervals 55479.801-55479.824 and 
55479.3652-55479.5483, respectively, [A], and during RXTE observa- 
tions covering the intervals MJD 55482.008-55482.043 [B] and MJD 
55483.715-55483.839 [C], respectively. The best-fitting models, defined as 
TBAbs* (nthcomp+gaussian) , are over-plotted as black solid lines. 
Points referring to XRT, JEM-X, ISGRI and PCA spectra are plotted as red, 
green, purple and blue points, respectively. In the lower panels, residuals 
with respect to the best-fitting models are plotted. 



tained. Since the total flux is F = Fs{l + e^), where Fs = 
aTg{Rs/d)^ is the observed flux coming from the colder region 
and y = 4(A;T'e/w-eC^)r^ is the Comptonization parameter, one 
has Rs = 1. 8(1 /kTsf^yF /(I + ev) dg.g km. Here, fcTg is ex- 
pressed in keV and ds.g is the distance to the source in units of 5.9 
kpc. By inserting the observed values (see Tab. |3) an emitting re- 
gion with radius between 3 and 15 km is estimated, indicating the 
NS surface as the most probable region of emission of the soft pho- 
tons. The iron line detected by PCA appears broad (cr ~1 keV) and 
intense (EW ~ 0. 2-0.3 keV), in accordance with the finding s of 
iMiUer et al.l ( l201lh and IChakrabortv & Bhattacharvval (120111) . Its 
energy is generally compatible with neutral iron (6.4 keV), even if 
larger values are also sometimes observed. 

The accuracy of the spectral modelling obtained from PCA 
observations was checked by considering the spectra obtained by 
XRT during the six Swift observations which partly overlapped or 
were very close in time to RXTE pointings. During all of these ob- 
servations XRT was operated in WT mode. Spectral fitting was per- 
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Figure 3. Error contours on the electron temperature and optical depth mea- 
sured during observations labeled as [A], [B] and [C] (see text and the cap- 
tion of Fig. [2), drawn considering the 90% confidence intei'vals calculated 
for two parameters (Ax^ = 4.61). 



formed by letting the relative normalisation of XRT with respect to 
PCA free to vary; values in the range 0.92-0.96, with a typical un- 
certainty of 0.01, were measured. The results obtained are given 
in Tab. |4] The combined spectra are well described by the model 
discussed previously, even if residuals are sometimes found at the 
overlap between the two energy bands. The best-fitting parameters 
of the simultaneous modelling of XRT and PCA spectra are plot- 
ted as red points in Fig. [T] and are visibly very close to the values 
obtained by modelling PCA spectra alone. The addition of XRT 
spectra then confirms that the spectral decomposition obtained by 
analysing PCA data alone is robust. The values of the absorption 
column obtained with the simultaneous modelling of XRT and PCA 
spectra are between 1.0 and 1.5x10^^ cm~^, compatible with the 
Chandra estimate. 

As it was noted earlier, the source emission is soft during 
the large part of the outburst. Nevertheless, the electron temper- 
ature measured during the first PCA observation (labelled as [B] 
in the top panel Fig. [T), 6.7(2) keV, suggests that a spectral tran- 
sition from a harder state took place during the outburst rise. To 
show this spectral change we plot in the top panel of Fig. |2] the 
unfolded spectra extracted duiing observations designated as [B] 
(covering the interval MJD 55482.008-5482.043) and [C] (MJD 
55483.715-55483.839; see top panel of Fig.[T). To further inves- 
tigate this spectral transition, we also considered observations per- 
formed by Swift and INTEGRAL before the start of the RXTE cover- 
age on MJD 55482.008 (see label [A] in the top panel of Fig.[T]l. We 
then fitted the 1-10 keV Swift XRT spectrum accumulated starting 
on MJD 55479.801, for an exposure of 2.0 ks (Obs. 00031841002), 
together with the 4.5-20 keV JEMX-7NTEGRAL and the 20- 
150 keV ISGRI-INTEGRAL spectra accumulated during Rev.975 
(spanning MJD 55479.365 - 55479.548, for an exposure of ~ 8.5 
ks). We normalised the flux observed by JEM-X and ISGRI to that 
observed by XRT by letting a normalisation constant free to vary. 
Given the relatively low count rates observed by JEM-X and IS- 
GRI, we held the normalisation constant between these two instru- 
ments fixed to I. No significant improvement in the modelling is 
obtained otherwise. The spectrum so produced is well fitted by the 
same Comptonization model used above, with best-fitting param- 
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Figure 4. Background-subtracted fractional amplitude of the first (top) and 
second (bottom) harmonic component used to model the pulse profile of 
J17480. Red triangles refer to profiles obtained by folding 0.2-10 keV XRT 
time series around i/q = 11.044885 Hz over the entire length of each 
observation, while black points are obtained by folding 300 s long segments 
of 2-60 keV PCA data around the same frequency. 



eters, kTs = 0.7(1) keV, kTe = 20tf keV and r = 4.4J:J;^ 
(see Table[5](. These indicate a hotter and optically thinner Comp- 
tonizing medium than what is observed during later observations. 
This is clearly shown in the top panel of Fig.|2l where the unfolded 
spectrum of the combined XRT, JEM-X and ISGRI spectra is la- 
belled as [A]. The best-fitting parameters of spectrum [A], as well 
as the luminosity extrapolated to the 0.05-150 keV energy band, 
are reported in Fig.[T|with green points. To better show the spectral 
transition performed by J17480 during the outburst rise, and to take 
into account the correlation between the values of electron temper- 
ature and optical depth measured with the Comptonization model 
we employed (see Eq.[T](, we plot in Fig.[3]the 90% confidence level 
error contours for these parameters. The contours were obtained by 
considering the surface for which the fit chi-squared was larger by 
Ax^ ~ 4.61 with respect to the best fitting value. 

A softening of the emission during the outburst rise is also 
indicated by the analysis of the spectra accumulated by JEM-X 
and by ISGRI on-board INTEGRAL, alone. We list in Tab. |6] the 
results obtained by fitting the spectra accumulated during each 
of the INTEGRAL revolution with the same model used so far, 
and keeping the relative normalisation between the two instru- 
ments fixed to one. Observations performed before MJD 55483 are 
well fitted by Comptonization in a medium with kTe > 10 keV, 
while subsequent observations indicate temperatures of the order 
of those indicated by the analysis of XRT and PCA spectra (~ 3 
keV). Such a transition is even clearer by taking into account the 
fluxes observed by JEM-X and ISGRI. Before MJD 55483 the ra- 
tio FjEMx (4.5 - 25.0) /FisGRi (20.0 - 150) is Ri 0.5 - 1, while 
it takes values in excess of 30 during observations performed when 
the source bolometric luminosity is larger Moreover ISGRI could 
not detect the source after MJD 55483, with typical upper limits 
of FisGRi{50 - 150) < few x 10""-10"^^ erg cm'^ s'^ 
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Figure 5. Background-subtracted profiles obtained by folding around 
VQ = 11.044885 the 2-60 keV PCA light curves accumulated during ob- 
servation labelled as [B] (spanning MJD 55482.008-55482.043), [C] (MJD 
55483.715-55483.839) and [D] (MJD 55487.565 - 55487.693), when the 
0.05-150 keV luminosity was 1.91(4) X lO^'^, 4.81(2) X lO^'^ and 
9.23(4) X 10 erg s~ , respectively. The difference of the phase of the 
various pulse profiles is interpreted in Sec. I3.3l and r3.4l in terms of the NS 
spin frequency derivative. Two cycles are plotted for clarity. 



(3 a confidence level), while the flux in the same energy band is 
~ 5 X 10"^" erg cm^^ s^^ during earlier observations (see the 
rightmost columns of Tab.|6}. 

Estimates of the bolometric luminosity were obtained by ex- 
trapolating the best-fitting models of PCA spectra to the 0.05- 
150 keV energy band, and using the mo st updated estimate of the 
distance to the cluster (d=5.9±0.5 kpc. iLanzoni et al.ll2010 '). The 
values we obtained are listed in Tab. [3] and plotted in panel (a) 
of Fig. [T] The maximum luminosity, Lmax ~ 9.23(4) x 10"^^ 
d| 9 erg s~^, is observed during the observation starting on i = 
MJD 55487.63(6) = 7.81(6) d since To, and labelled as [D] in the 
top panel of Fig.[T] To give an estimate of the X-ray flux before the 
start of the RXTE coverage of the outburst, we consider XRT and 
INTEGRAL observations previously labelled as [A], yielding an 
estimate Lmin = 0.44(6) x 10^'' dig erg s"\ 



3.2 The pulse profile 

To analyse the properties of the pulse profile of J17480 we consid- 
ered the time-series recorded by PCA in the 2-60 keV energy band 
and by the XRT in the 0.2-10 keV band. We first applied barycen- 
tric corrections on the times of arrival of photons recorded by PCA 
and XRT, to take into account the motion of the spacecrafts in the 
Solar System. To this end we considered the position, RA=17^M8™ 
4': 83 1( 4), DEC=-24° 46' 48" 87(6), determined by [Poolev et al.l 
HQIu) by comparing a Chandra observation of the source during 
the outburst with a deeper past Chandra observation of Ter zan 5, 
performed when J17480 was quiescent (Heinke et al. 2006). This 
postion is compatible with the one derived by Riggio et al. (2013) 
from a Moon occultation of the transient, observed by RXTE dur- 
ing the first days of the outburst. The NS motion in the binary sys- 
tem was corrected by using the orbital parameters found by PI 1, 
who performed a timing analysis on a subset of the observations 
considered here. 

PCA time series were split in 300 s long intervals and folded 
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Figure 6. Energy dependence of the RMS amplitude of the pulse pro- 
file shown by J 17480 during PC A observations performed during the in- 
terval MJD 55482.008-55482.043 (labelled as [B]) and MJD 55483.715- 
55483.839 (labelled as [C]), respectively. 



in 16 phase bins aro und vp = 11. 044885 Hz (see Pll). We adopt 
the criterion given bv lLeahvl lll987h to determine a detection thresh- 
old. A signal was detected with a significance larger than 3 cr in all 
PC A observations. Given its lower effective area, the count-rates 
recorded by XRT are significantly lower and the whole length of 
each observation was considered in the folding procedure to max- 
imise the counting statistics. A signal could be detected only in 6 
out of 18 XRT observations; however, when a detection was miss- 
ing the upper limit on the signal amplitude was larger than the am- 
plitude indicated by the analysis of quasi-simultaneous PCA obser- 
vations. The limited signal-to-noise ratio of XRT data sets is then 
the most plausible reason for these non-detections. 

A two-components harmonic decomposition well describes 
the observed pulse profiles. The evolution of the background- 
subtracted 2-60 keV fractional amplitude of the first {Ai) and sec- 
ond harmonic (^2) observed by PCA is plotted as black points in 
the top a nd bottom panel of Fig. ID respectively. As it was already 
noted by Papit to et al.l J201lh . a sudden variation in the pulse pro- 
file properties takes place during the outburst rise. Most noticeably, 
the fractional amplitude of the first harmonic observed by PCA de- 
creases from values of ~ 27%, measured before MJD 55483, to 
~ 2% during later observations. The second harmonic amplitude 
is also observed to decrease, though much more smoothly than 
the amplitude of the fundamental. During the decaying stage of 
the outburst, the amplitudes of the two harmonic components do 
not show any other step-like behaviour, following a smooth trend. 
Even though the different instrumental responses do not allow an 
immediate comparison, also the amplitude observed by XRT in the 
0.2-10 keV energy band follows a similar trend (see red triangles 
inFig.|4ll. 

To show the evolution of the pulse profile shape taking place 
around MJD 55483 we plot in Fig. [5] the background-subtracted 
profiles observed by PCA during the observation labelled as [B], 
[C] and [D], covering the intervals MJD 55482.008-55482.043, 
MJD 55483.715-55483.839 and MJD 55487.565-55487.693, and 
during which the 0.05-150 keV X-ray luminosity was 1.91(4) x 
10^^, 4.81(2) X 10^^ and 9.23(4) x 10^^ erg s"\ respectively. 

Also the energy dependence of the pulse amplitude changes 
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Figure 7. Background subtracted pulse profiles in four energy bands accu- 
mulated during PCA observations of group [B] (see text). The solid curves 
are the respective best-fitting harmonic decompositions. Two cycles are 
plotted for clarity. 



during the outburst rise. To show this, we consider the observations 
performed by RXTE during the intervals [B] and [C]. Such a choice 
is made since the pulse properties observed during group [C] data 
qualitatively represent what is observed during the rest of the out- 
burst. The energy dependence of the (background subtracted) pulse 
RMS amplitude, r{E) = [{Aj + A\)l2\^l'\ is plotted in Fig.|6] 
During observations labelled as [B] the pulse RMS amplitude is 
roughly constant with an average value of 0.19. After MJD 55483 
instead, the pulses are nearly completely suppressed at low ener- 
gies and become stronger as the energy increases up to 30 keV (see 
curve labelled as [C] in Fig.|6). While during observations of group 
[B] the profile is nearly sinusoidal and shows the same shape in 
different energy bands, the phase of the two harmonics becomes 
slightly energy dependent after MJD 55483. To see this, we plot in 
Fig.|7]and[8]the background subtracted pulse profiles in four PCA 
energy bands, accumulated during observations labelled as [B] and 
[C], respectively. During observations of group [B] the profile has 
a small harmonic content and shows a similar shape in different en- 
ergy bands. During observations [C] the phase difference between 
the first and the second harmonic is negligible at low energies, re- 
sulting in a nearly symmetric profile shape, whereas, at higher ener- 
gies, the fundamental lags the second harmonic resulting in a more 
skewed profile. 



3.3 Timing analysis. 

The spin evolution of the source during the outburst was investi- 
gated under the assumption that the evolution of the phase of the 
pulse profile is a good proxy of the NS spin, that is 
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Figure 8. Background subtracted pulse proiiles in four energy bands accu- 
mulated during PCA observations of group [C] (see text). The solid curves 
are the respective best-fitting harmonic decompositions. Two cycles are 
plotted for clarity. 



-^(t) - -/-o = 



Wit') - vf]dt', 



(2) 



To 



where p(t) is the spin frequency of the NS, Vf is the fre- 
quency around which the time series are folded and To = MJD 
55479.81229 is the reference epoch of the timing solution. The 
phases calculated over the first and the second harmonic compo- 
nent are plotted in the top panels of Fig. |9] and [TO] respectively. 
Phases computed over the profiles observed by XRT in the 0.2-10 
keV do not show a signficant scatter with respect to PCA phases 
(2-60 keV). This can be seen from the bottom panels of Fig.|9]and 
1101 where we plotted the residuals of XRT (red triangles) and PCA 
(black points) phases with respect to the best-fitting models defined 
below. We conclude that the influence on the observed phases of the 
different energy band and response of the two instruments is small, 
and consider both the values measured by XRT and PCA. Con- 
sidering PCA phases alone would not significantly alter the results 
obtained. 

To characterise phenomenologically their evolution we fit 
them with the relation 



4>(t) - 00 = PN{t - To) + Rorb(t - To). 



(3) 



Here, Pjv(f — To) = X/n=o ^nPn(t — To) is a A''-th order polyno- 
mial expansion, and the p„'s are discrete orthogonal polynomials 
over the considered data set. Such a choice makes the various fitted 
parameters, c„, independent on the d egree of the fitted polynomial, 
A*' (see, e.g.. lBevington & Robinsoij i2003). Defining Xi = ti — To, 
the terms of the polynomial expansion are obtained through the re- 
currence relation, p„+i(a;) = {x — A„+i)p„{x) — B„p„-i{x), 
withpo(2:) ~ 1 andp_i(a;) — 0. The set of constants [An, Bn) is 
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Figure 9. Phases computed over the first harmonic component by folding 
PCA (2-60 keV; black points) and XRT (0.2-10 keV; red triangles) time 
series around vo = 11.044885 Hz, together with the best fitting 8-th order 
polynomial (see text) plotted as a solid line (top panel). Residuals in units of 
sigma with respect to the best fitting model are plotted in the bottom panel. 



determined by imposing orthogonality between the various terms 
of the polynomial expansion, i.e. "Y^ .{1 / a^)pk{xi)pj{xi) = 0, for 
k 7^ j, k,j = 0, ..., N. The weights 1/af simply reflect the uncer- 
tainties on the various phase points. The term Rorbit) is added to 
Eq.[3]to take into account the effect on the phases of a small dif- 
ference between the orbital parameters used to correct the time se- 
ries (the orbital period Porb, the projected semi-major axis a sini 
and the epoch of zero mean longitude T*) and the actual orbital 
parameters of the system. An expression for Rorbif) is given by 
lDeeteretal1 ( ll98lh . 

Small corrections to the orbital parameters are found, lead- 
ing to the estimates given in Table. [T] (obtained considering a fit 
with A'' = 8, see below). The relatively large values the Xr ob- 
tained from the phase modelling indicates the presence of phase 
noise. Slight differences are found between the orbital parameters 
obtained by fitting the two different harmonics, indicating that such 
a noise induces phase variations also on a timescale comparable 
with the orbital period of the system. To give conservative estimates 
of the orbital parameters we then quote the average among the pa- 
rameters obtained from the phases of the first and second harmonic. 
The uncertainty on each parameter is estimated to encompass the 
1 a confidence level interval of the value indicated by the analysis 
of each component. The values of the orbital parameters we obtain 
are compatible with the estimates obtained by Pll on a reduced 
data set. The new set of parameters is then used to correct the time 
series. The errors induced on the phase estimates by the uncertain- 
ties in the orbital parameters [see Eq. 3 in Burderi et al. (2007)] are 
now summed in quadrature to the statistical phase uncertainties. By 
doing so any residual orbital modulation of the pulse phases due to 
the uncertainty affecting the orbital parameters is already included 
in the phase errors. The fitting of the phases is then repeated by 
using Eq.|3]with Rorbit) = 0. This makes the function defined by 
Eq. [3] eventually orthogonal on the considered data set. 
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Figure 10. Same as Fig.|9]regarding the phases computed over the second 
harmonic component. 



Starting from A^ = 1, the order of Pjv(t — To) is increased 
and the phases re-fitted with the higher-order polynomial until the 
improvement in the x^ of the fit is not significant on a 3 ct con- 
fidence level according to an F-test, after two consecutive order- 
increments. Despite that such a criterion is somewhat arbitrary, the 
orthogonality between the various orders of the polynomial expan- 
sion grants that the estimates of the fitted parameters, Cj, j ^ A*', 
do not depend on the polynomial order, allowing a coherent phe- 
nomenological description of the phase evolution. 

A by-eye inspection of Fig. |9] and Fig. 1 101 already shows the 
need for a quadratic term to fit the phase evolution. This is further 
confirmed by the extremely large values of the F function (> 10^) 
obtained when adding such a term to a simple linear function. Un- 
der the assumption that the phases are a good tracker of the NS spin 
(Eq.|2](, we have 



mt) 






(4) 



dt'-"') 

which shows how the best-fitting value of every parameter c„ trans- 
lates into an estimate of the average value of the (n-l)-th derivative 
of the spin frequency. The values of the spin frequency at To and 
of the average spin frequency derivative we obtain from the fit of 
the first and the second harmonic phases are quoted in Table |2] 
The two estimates of the average spin frequency derivative are very 
close to each other, and a value of < !> >= 1.48(2) x 10~^^ 
Hz s~^ encompasses both. By considering a smaller dataset (the 
first 25 days of the outburst compared to the 40 days considered 
here), Pll obtained values of the average spin-up rate of 1.22(1) 
and 1.68(l)xl0^^^ Hz s^^, from the phases computed over the 
first and the second haiTnonic, respectively. This clearly shows how 
the considered dataset may influence the determination of the spin 
parameters of an accreting pulsar when profile it shows is affected 
by timing noise. However, the scatter between the values obtained 
from the first and second harmonic greatly decreases when a longer 
interval is considered, suggesting how the estimate given here is 
more reliable. 



Figure 11. Black solid lines represent the second derivative of the 8-th or- 
der polynomial that best-fits the phases computed over the first (top panel) 
and the second harmonic (bottom panel). Dashed black curves mark the re- 
spective 1 cr confidence-level intervals. The change of sign of the second 
derivative of the best-fit polynomial is done to make the plotted function 
equal to the spin frequency derivative (see Eq.|4) and to facilitate compari- 
son with light curve plotted in the top panel of Fig.[T]and reproduced in the 
bottom panel of this figure by using red circles (flux units are arbitrary). 



The phase evolution is anyway far from being adequately 
fitted by a simple quadratic function and significant improve- 
ments of the description are successively obtained adding higher 
order terms, up to A'^ — 8. Considering that the phases com- 
puted over the two harmonic components have similar uncertain- 
ties, the quality of the modelling of the second harmonic phases 
(xl = 1371.8/673 = 2.04) is significantly better than for the 
first one (x? = 4614.2/760 = 6.07). Yet, these values indi- 
cate the presence of phase fluctuations which are not described 
by a polynomial expansion of a relatively low order. These fluc- 
tuations affect in a different manner the two harmonic compo- 
nents, though, and the trend followed by second harmonic phases 
is smoother than that followed by the phases computed over the 
first harmonic. To quantify this we evaluated the total RMS am- 
plitude of the phase residuals with respect to the best-fitting poly- 
nomial as the squared sum of the statistical uncertainty on the 
phase measurement (which depends on the amplitude of the con- 
sidered harmonic and on the count rate) and of the intrinsic phase 
noise, cr^ ,,„, — cyk,,„tat^l^,.noise^ where k is the harmonic num- 
ber (see iHartman et al.ll2008al and references therein). By taking 
an average value for the statistical uncertainties one obtains that 

0-l,noise = 0.055 and 0-2, noise — 0.017. 

An enhanced stability of the second harmonic phases is also 
indicated by the smooth trend followed by the second derivative 
of the best-fit polynomial, which equals in modulus the spin fre- 
quency derivative under the assumption that the NS spin is well 
tracked by the pulse phases (see Eq.|4j. In the top and bottom panel 
of Fig. \TT\ we use black solid lines to plot the second derivative 
of the best-fitting function of the phases computed over the first 
and second harmonic, respectively. The curve computed over the 
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Table 1. Orbital parameters of J17480 obtained by modelling with Eq.[3] 
the phases computed on the first (left column) and on the second (central 
column) harmonic, together with the x^ of the respective modelling. Uncer- 
tainties are given at 1 cr confidence level and were scaled by a factor yJXr 
to take into account that x^ > 1- The average between the values obtained 
consideiing the two harmonic components are listed in the rightmost col- 
umn, with the uncertainties estimated conservatively such as to overlap the 
values indicated by the two components. 





P* harmonic 


2"'' harmonic 


Average 


a sin i/c (It-s) 


2.4972(2) 


2.4972(1) 


2.4972(2) 


Porb (s) 


76588.4(1) 


76587.92(6) 


76588.1(3) 


T* (MJD) 


55481.78038(2) 


55481.78052(1) 


55481.78045(8) 


e 


< 6 X 10-'' 


< 1 X 10-* 


< 2 X 10-4 



5192.99/757 



Table 2. Frequency at the reference epoch and average spin frequency 
derivative obtained by fitting the first and second hamionic phases with the 
8-th order polynomial defined by Eq.|3](with Rorb{^ — Tq) = 0), together 
with the relative x^ and amplitudes of residuals (see text). Uncertainties are 
given at Icr confidence level and were scaled by ^/xr to take into account 
that Xr > 1- 



1436.74/668 



!/(ro)(Hz) 

<;>>(Hzs-i) 


V* harmonic 

11.0448803(7) 
1. 489(4) X 10- 12 


2"'' harmonic 

11.0448797(5) 

1.468(2)xl0-i2 


'^rms 
(Jstat 
^noise 


0.059 
0.022 
0.055 


0.024 
0.017 
0.017 


x^ 


4614.2/760 


1371.8/673 



first harmonic phases shows an oscillatory behaviour, especially to- 
wards the boundaries of the available data set. Such rapid variations 
of the spin up rate in absence of simultaneous swings of the X-ray 
light curve (the shape of which is reproduced in the bottom panel 
by using red circles; see also the top panel of Fig. [T) indicate how 
hardly the phase variations computed over the fundamental fre- 
quency component can be interpreted in terms of a variable torque 
acting onto the NS (see Sec. |4.3t . This possibly indicates that the 
amplitude of the timing noise affecting the phases evaluated over 
this component is larger than the deviations from a parabolic trend 
due to the non constancy of the accretion torque, making the lat- 
ter practically unmeasurable from this harmonic component. The 
phases computed over the second harmonic have instead a much 
smoother behaviour. The spin up rate increases during the first days 
reaching a maximum value of2.8(l)xl0~i2 Hz s-^ at the epoch 
6.0±0.5 d since To, and subsequently smoothly decreases until it 
becomes ill-defined towards the boundary of the given data-set. 
Comparing the spin up rate calculated over the second harmonic 
phases to the X-ray light curve, it appears clear that the vii) curve 
follows much more closely the expected dependence on the instan- 
taneous mass accretion rate. 

These results indicate how the phases computed over the sec- 
ond harmonic are possiby a better tracker of the spin frequency of 
the NS. For this reason we use only this component to investigate 
if the observed behaviour is compatible with the accretion of Kep- 
lerian angular momentum of the matter at the inner edge of the disc 
(see next section). This choice is further motivated by the analo- 
gous behaviour observed in other accreting pulsars, where the sec- 
ond harmonic phases are often observed to be less affected by tim- 
ing noise than the fundamental phases (see the Discussion, where a 
plausible explanation for such a behaviour is also proposed). 



3.4 Torque modelling 

To analyse whether the evolution of the pulse phases computed over 
the second harmonic can be described in terms of the torque im- 
parted to the NS by the accreting matter, we assume that the torque 
depends on a power, /3, of the mass accretion rate, 27r/, i/ oc M^ , 
and that the mass accretion rate is well tracked by the X-ray lumi- 
nosity, M — RtLx/riGM,, where A/,, R, and /, are the mass, 
radius and moment of inertia of the NS, respectively, and 77 ~ 1 is 
the accretion efficiency. We thus consider a relation 



u{t) = u{i) 



Lx{t) 
Lx{i) 



(5) 



to express the dependence of the spin frequency derivative on the 
accretion luminosity. The bolometric accretion luminosity is es- 
timated by extrapolating to the 0.05-150 keV energy band the 
best-fitting spectral models discussed in Sec. 13.11 The epoch t — 
MJD 55487.63(6) = 7.81(6) d since To, is the epoch at which the 
maximum luminosity, Lxit) ~ Lmax = 9.23(4) x 10''^ d| 9 erg 
s-i is attained. 

If the folding frequency is close to the actual spin frequency 
of the source and the frequency variation throughout the outburst is 
small, Eq.|2]can be expressed as 

0(t) - 00 ==; -5,y {t - To) - /^^ dt' J^^ v{t") dt" = 

= -Si^ {t - To) - !>(i) /^^ dt' J^^ [Lx{t")/Lx{t)f dt",(6) 

where Su — J^(To) — J^/, and we have inserted the assumed depen- 
dency of the frequency derivative on X-ray luminosity (Eq.|5j. One 
thus obtains a relation that can be used to fit the observed phase 
evolution depending on the parameters (po, i^{Tq), iy{t) and /?. 

In practise, we approximated the luminosity profile of the out- 
burst Lxit) by using a natural smoothing cubic spline, which is 
plotted as a solid curve in panel (a) of Fig. [T] We then performed 
numerically the double integral of the function L^{t) for a set 
of values of /3. We thus obtained a set of functions that were 
used to fit the observed evolution of the phases computed over 
the second harmonic. The best-fit is found for /3 — 1.07 ± 0.03 
with a x^ = 1574/672 = 2.34, a value not far from the x^ 
of the best-fitting 8-th order polynomial (Xr ~ 2.04). This sug- 
gests a link between the phase variation computed over the sec- 
ond harmonic and the X-ray flux measured from the source. The 
best-fit values of the spin frequency at the beginning of the out- 
burst To and the spin frequency derivative at the outburst peak are, 
iy(To) = 11.044884819(6) Hz and i>{i) = 2.755(5) x 10"^^ Hz 
s"'^, respectively. The errors on the fitted parameters (quoted at the 
la confidence level) were evaluated by scaling the statistical errors 
by yXr to take into account the fact that Xr > 1- The best-fitting 
model is plotted over measured data points, together with residuals, 
in Fig.[l2] 

To interpret the timi ng noise affec ting the pulse phases ob- 
served from AMSP, Patrun o et al .1 ( 120091) proposed that a large frac- 
tion of the unmodelled phase variability could be ascribed to a cor- 
relation between the phases and the X-ray flux. This would possi- 
bly reflect a dependence of the hot-spot location on the accretion 
rate. To check this hypothesis in the case of J 17480 we substi- 
tuted the double integral of the X-ray luminosity in Eq. |6]with a 
linear term in the X-ray luminosity, k{t) — D Lxit). The value 
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described by Comptonization in a hotter and more optically trans- 
parent medium (with an el ectron temperat ure of ~20-60 keV and 
optical depth ~0.7-2.5, see |Poutanenl200q . and references therein), 
not showing strong spectral variability during an outbursts. Only 
during the earliest observations, when the source was fainter than 
~ 3 X 10''^ d§ g erg s^^, the spectrum observed from J17480 some- 
what resembles those displayed by AMSPs. At larger luminosities 
instead, the observed spectrum is much more similar to those show n 
by brighter, non-pulsing, accretors (see, e.g. . Barret & Olivel2002h . 
The decrease of the electron temperature of the Comptonized spec- 
trum as the accretion luminosity increases can be interpreted in 
teiTns of an increase of the flux of soft photons which efficiently 
cool the Comptonizing region. In this context, it is worthwhile to 
note that based on the anal ysis of its c orrelated colour and a periodic 
timing variability, Altamir ano et al.l (12010) and Chakrabor tv et al.l 
( l201l[) concluded that this source transited from an Atoll-like to a 
Z-source like behaviour on MJD 55485, when the X-ray luminosity 
was close to attain its peak value. The hard to soft spectral evolution 
reported here is then part of a more general change of the properties 
of the accretion flow near the NS surface as the mass accretion rate 
increases. 



Figure 12. Best fitting model according to the physically-motivated torque 
model defined by Eq.[5]and|6]of the second harmonic phases (top panel), 
together with residuals in units of a (bottom panel). 



of the chi-squared we obtain by fitting this model, Xr ~ 35.8, 
clearly indicates that the observed phase variability cannot be ex- 
plained in terms of such a correlation alone. Nevertheless a phase- 
flux correlation may still explain a small fraction of the phase 
variability, in addition to the accretion rate dependent torque. By 
adding k[t) to Eq.[6]the chi-squared of the fit decreases in fact by 
Ax^ = 70, with one degree of freedom less. The best-fitting pa- 
rameters we obtain after the addition of this linear term to the model 
are, v{Tq) = 11.044884842(6) Hz, v{t) = 2.602(6) x 10"^^ Hz 
s-\ P = 0.97 ± 0.04 and D = 0.007(1) cycles s erg"\ 



4 DISCUSSION 

4.1 The energy spectrum 

In Sec 13.11 we presented an analysis of the spectra accumulated 
by PCA on-board RXTE (2.5-30 keV), by XRT on-board SWIFT 
(1-10 keV) and by the JEM-X (4.5-25 keV) and ISGRI (20-150 
keV) detectors on-board INTEGRAL. The phase-averaged contin- 
uum spectrum of J 17480 is well described by thermal Comptoniza- 
tion. During the first few days of the outburst the temperature of 
the Comptonizing medium is kT^ ^ 10 keV. It then decreases to 
values ~ 3-3.5 keV as the source bolometric luminosity exceeds 
^ 3 X 10^^ ds 9 erg s"\ on MJD ~ 55483, without showing 
further strong deviations from these values. Correspondingly, the 
optical depth increases from r ~ 5 to 10. Unfortunately, because 
there is no observational coverage of the later stages of the out- 
burst, it is not possible to ascertain the presence of a correspond- 
ing soft-to-hard transition and at which luminosity does it happen. 
The spectral transition observed during the outburst rise, as well as 
the softness of the spectrum characterising the emission of J17480 
during most of the observations, is not typical of faster accreting 
pulsars. If compared with J17480, AMSP attain lower peak lumi- 
nosities ( ^ few xlO'^'' erg s~^), and their spectrum is generally 



4.2 The pulse profile shape 



We presented in Sec. |3.2| an analysis of the 1 1 Hz pulses shown by 
J17480, based on PCA and XRT observations. The pulse shape can 
be modelled by using no more than two harmonic components. A 
major change of the pulse properties is observed at MJD ^ 55483, 
while the source flux is increasing and the spectrum softening; the 
fractional amplitude of the first harmonic component suddenly de- 
creases from ~ 27% to ~ 2%, attaining a value of the same order 
of that of the second harmonic (see Fig.Q. The drop in the pulse 
strength is higher at soft energies, since the pulse amplitude be- 
comes strongly energy dependent, showing a marked increase of 
the amplitude with the energy of the photons up to 30 keV (see 
Fig.|6]and[8l(. 

The shape of the pulses produced from the surface of an ac- 
creting rotating NS was modelled by a number of authors (see, e.g., 
Muno et al.l |2002|: [Poutanen & Beloborodovl l2006l : ICadeau et all 
20071 : iLamb et al.l 12009! . and references therein) and depends on 
the geometry of the system (i.e., on the relative angles between 
the spin axis, the magnetic axis and the line of sight), on structural 
properties of the NS (mass, radius, shape and spin), on the beaming 
pattern of the emitted radiation, as well as on the size, temperature 
and brightness contrast of the accretion cap with respect to the rest 
of the NS surface. During most of the outburst the pulses shown 
by J 17480 are of low amplitude (few per cent) a nd simple, simi- 
lar to the majority of AMSP (see lLamb et al.ll2009!. L09 hereafter). 
However, the amplitude decrease by more than a factor of 10 tak- 
ing place on MJD ~ 55483 on a timescale of less than one day 
(see Fig. |4j, indicates how some of the properties that determine 
the pulse amplitude must have changed significantly as the X-ray 
luminosity was increasing during the outburst rise. Variations of the 
pulse amplitude of a similar amount were already observed from a 
number of AMXPs; swings f rom 3 to 25 per cent have been ob- 
served from XTE J1807-294 dZhang e t al. 2 006), and front 4 to 1 4 
per cent in the case of SAX J1808.4-3658 OHartman et al.ll2008bh . 
In those cases, though, the amplitude variations were closely re- 
lated to flaring episodes during the outbursts. On the other hand, 
the decrease of the pulsed fraction shown by J 17480 is step-like 
and no significant correlation holds between the pulsed fraction and 
the X-ray flux afterwards. In the following a number of hypotheses 
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to account for the observed variation of the pulse properties during 
the outburst rise is discussed. Unfortunately, a coverage of the later 
stages of the outburst decays is not available, preventing to investi- 
gate if a similar transition happens as the source dims and at which 
luminosity would it take place. 

A change of the spot size is hardly the main driver of the ob- 
served pulse fraction decrease; the spot size is in fact expected not 
to vary by more th an ^ 10% as the luminosity doubles, and the 
results obtained by i Muno et al.l 1 120021) and L09 show how the de- 
pendence of the pulse amplitude on the spot size is much weaker 
than it would be needed to explain the observed variation. 

Pulse amplitude variations can be explained in terms of move- 
ments of the hot spot location on the NS surface, possibly caused 
by the in-falling matter getting attached to a different set of mag- 
netic field lines. An increase of the spot latitude can be particularly 
efficient in decreasing the pulse amplitude, especially if the angle 
between the magnetic and the spin axis is small. This was invoked 
by L09 to explain some of the properties of the pulses shown by 
AMSPs, such as small and variable amplitudes, and unstable pulse 
phases. According to their results a change of roughly 15 — 30° 
of the spot latitude towards the magnetic pole would be able to 
produce an amplitude decrease of the order of that observed from 
J 17480. Further, under the hypothesis that the two spots on the NS 
surface are antipodal, a re-sizing, or shift, of the impact region of 
the in-falling matter on the NS could increase the fraction of rota- 
tional phases during which both spots on the NS surface become 
visible. This would indeed contribute to explain the observed re- 
duction of the pulse amplitude and strengthening of the second har- 
monic amplitude relative to the first one. Although this scenario is 
then in principle able to explain the decrease of the pulse amplitude 
integrated over the entire PCA energy-band, it gives no sufficient 
clues for the simultaneous change of the energy dependence of the 
pulse profile shape. 

An amplitude variation of the order of that observed could 
also be produced by a significant change of the angular distri- 
bution of the photons that reach the observer, possibly triggered 
by th e source luminosity crossing a certain threshold level (see, 
e.g., IParmar et al. 1989, who interpreted as such the luminosity 
dependent changes of the pulse profile of the 42 s high-field pul- 
sar, EXO 2030+375). While the radiation transmitted through a 
slab is strongly peaked along the normal to the NS surface, scat- 
tered radiation is expected to peak at intermediate angle s, in a fan- 
like emission pattern (e.g. Sunvaev &Titarchuklll985l) . The rela- 
tions given by IPoutanen & Beloborodovl ( 120061 see, e.g., Eq. 51 
and 61) show that an amplitude change of the same order of that 
observed can be reproduced if the angular distribution of the emit- 
ted photons switches from being strongly peaked along the normal 
[I (a) oc (1 + hcosa), where a is the angle to the normal and 
h ~ 2], to a fan like distribution (h~ -1). However, such a dras- 
tic change of the pattern of the emitted radiation would produce a 
change of ~ 0.5 in the pulse phase, which is not observed. 

The decrease of the pulse amplitude shown by J 17480, as well 
as its spectral hardening, take place roughly simultaneously to the 
spectral softening of the phase-averaged emission discussed in the 
previous section. It is then intriguing to explain these transitions 
in terms of a variation of the properties of the matter flow close 
to the NS surface. The observed behaviour would be in fact ex- 
plained if a large fraction of the source luminosity, and in partic- 
ular its softer part, becomes not pulsed during the outburst rise. 
This could happen because the inner rim of the disc approaches 
the NS surface, with an increased part of the disc non-pulsed emis- 
sion falling in the considered energy band, and/or because an in- 



creasing fraction of the in-falling matter is not channelled to the 
magnetic poles, but rather accretes more evenly at the NS equator 
in a sort of optically thick boundary layer. The properties of the 
phase- averaged spectrum presented in Sec. 13. II would favour the 
latter hypothesis. The possibility that not all the plasma from the 
accretion disc is lifted off by magnetic field lines at the magne- 
tospheric boundary, with a fraction penetrating between the mag- 
netic field lines via a Rayleigh-Taylor instability or because of in- 
complet e coupling of the plas ma a nd the field l ines w as studied 
by, e.g., ISpruit & Taiuiil ( Il990l) and iMiller et al.l ( 1 19981 see refer- 
ences therein, also). Recent 3D MHD si mulations of disc accre- 
tion onto a magnetised NS performed by iRomanova et al.l ( l2008h 
showed how a transition to this twofold accretion channel may 
take place when the mass accretion rate becomes larger than ~ 
0.1{B/W^GfMEdd- Above this level it could then be that only 
a fraction ~ 10% of the accreted mass is channelled to the poles. 
The rest would be accreted preferentially at the star equatorial belt 
producing the observed optically thick Comptonized component, 
reminiscent of emission from non pulsing NS in LMXB. The hard- 
ening of the energy dependence of pulse amplitude would be nat- 
urally explained by this scenario. The pulsed emission originating 
from the polar caps is in fact expected to be harder than the non- 
pulsed emission produced in a denser equatorial boundary layer. 
Interestingly, an increase of the pulse amplitude with energy simi- 
lar to that of J 1 7480 was observed o nly from th ree AMSPs, SWIFT 
J1756.9-2508 jPatruno et al]|2010l), Aal X-1 JCasella et al.ll2008l) 
and SAX JI748.9-2021 jPatruno et aLll2009l) . the last two being 
sources that showed pulsations only intermittently. To conclude, it 
has to be noted that a pulse amplitud e increase at ene rgies Si 50 
keV was also observed by Falanga & TitarchukI (120070 from IGR 
J00291+5934. Since it concerned high photon energies, the authors 
interpreted such an increase as due to the decrease of the electron 
scattering cross section with energy. This effect clearly does not ex- 
plain the pulse amplitude properties observed from J 17480, as well 
as from the other mentioned sources, since the pulse amplitude is 
observed to increase even at low ( ^ 2.5 keV) energies. 

4.3 Phase evolution. 



We presented in Sec. |3.3| and |3.4| an analysis of the temporal evo- 
lution of the phase of the pulse shown by J 17480 during its 2010 
outburst, based on the whole data set available of PCA and XRT ob- 
servations. The results obtained confirm the earlier finding that the 
NS in this system accelerates while it accretes mass (Pll). While 
the analysis performed on the phases computed over the first and 
second harmonic gives slightly different values of the average fre- 
quency derivative, a value of < !> >= 1.48(2) x 10~^^ Hz s~^ 
encompasses both estimates. 

The observed frequency derivative of the signal emitted by 
a pulsar is the sum of the intrinsic change of the NS spin and of 
the acceleration of the pulsar with respect to the observer. While 
the apparent spin frequency derivative due to ce ntrifugal acceler- 
ation ( < ly^'^d/c ~ 1.4 x 10"^* ds.g Hz s~ MShklovskiilll970l. 
where /i = 85 ± 10 km s~^ is the velocity of the proper motion 
of Terzan 5 estimated bv lFerraro et al .1120091) and to galactic accel- 
eration (~ 2 X 10^^^ Hz s~^, obtained by ins erting the galactic 
coord inates of J 17480 in the relation given by lOamour & Tavlon 
1 19911) can be safely neglected, the acceleration of the pulsar in the 
gravitational fie ld of the cluster can in principle have a larger effect. 
IPhinnevI < 1 9931) estimated the maximum acceleration along the line 
of sight for an object located at a projected angular distance 6± 
from the cluster centre as l.lu GMcyi{< 9±) / (ctv d^ O^^) Hz s^^ , 
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where Mcyi{< 0±) is the mass within a projected distance smaller 
than d 9± from t he cluster centre. Con sidering the cluster parame- 
ters estimated bv lLanzoni et alj ( 120101) from the surface brightness 
profile of Terzan 5 (core and tidal radius of 9" and 277", respec- 
tively, and a total bolometric luminosity of 8 x 10'' Lq) and the 
projected distance of the pulsar from the cluster centre (6± ~ 4") 
to integrate the King density profile ( King 1962 ) up to 6±, yields 
a mass Mcyi(< 0±) ~ 5 x 10* M0, where a mass-to-light ra- 
tio of 3 was assumed. Inserting such a value in the relation given 
above yields an apparent frequency derivative of ~ 6 x 10~^^ Hz 
s~^, which is more than two orders of magnitude lower than the 
observed derivative and can be safely neglected. 

The average derivative of the signal frequency is compatible 
with the average acceleration imparted to the NS by the accretion 
of the Keplerian angular momentum of disc matter. By considering 
only this torque, A'^o, and assuming that the X-ray luminosity tracks 
the mass accretion rate (see Sec. l3.4t . the expected spin frequency 
derivative of an accreting NS is 



_No_ 

2TTh 



R, 



2-Kh 



2n I,r]VGAi: 



(7) 



where Ri„ is the inner disc radius. By inserting the measured value 
of < !> > in Eq. [T] an estimate of the average inner disc size of 
< R,„ >~ 64 7?2 7?;jo Af.,1.4 /|5(< L:„ > / 5)-" km is 
obtained, where i?«,io is the NS radius in units of 10 km, M^.ia 
the NS mass in units of 1.4 Mq, It, 4,5 the NS moment of inertia in 
units of 10*^ g cm^ and L37 is the accretion luminosity in units of 
10"^^ erg s~^ . Such a value is clearly compatible with the upper and 
lower bounds set by the presence of pulsations, 7?« < Rin ^ Re, 
where Re = 338 A/^ \^ 4 km is the corotation radius around the 
NSinJ17480(Pll). 

The analysis presented in Sec. 13.31 shows that a constant fre- 
quency model is not compatible with the observed evolution of the 
phases of J 17480. Yet, values of the reduced x^ in excess of one 
are obtained even if an 8-th order polynomial expansion is consid- 
ered as a model, indicating the presence of timing noise. Before 
discussing the results obtained by modelling the observed phase 
evolution with a physically-motivated torque model, the influence 
of such a noise on the phase determination has to be discussed. 

Noisy variability in the phases of accreting pulsars is not new 
and manifests itself through the presence of unmodelled phase 
residuals with respect to the considered timing solution. In partic- 
ular, it affects the pulses shown by AMSPs a nd it was subject of 

3 



a profound analysis and discussion (see, e. g. . iBurderi et alj|200€ ; 



2008J ; iRiggio et alj |2008£ : 



Papitto et al.' '2007*; 'Hartman et al.' 
Lamb et al. 2009; Patruno et al. 2009i, and references therein). The 
analysis of the temporal evolution of the phases of the pulse pro- 
file shown by a pulsar allows in principle the most accurate possible 
knowledge of the rotational evolution of the NS. However, irregular 
variability of the pulse phases not related to intrinsic spin changes 
of the NS hinders this possibility. Indeed, this is particularly trou- 
blesome for transiently accreting pulsars such as AMSP which are 
expected to experience relatively small frequency variations during 
their outbursts. In those cases the amplitude of the phase variability 
induced by timing noise can be comparable to that produced by a 
spin frequency derivative. This is certainly not the case of J 17480, 
since, despite the presence of timing noise, the pulse phase is ob- 
served to perform more than one phase cycle across the outburst 
(see Fig.[9landll0t. This source can be therefore regarded as a par- 
ticularly relevant case to the study of accreting pulsars, since the 
the contributions of timing noise and of the rotational evolution to 
the pulse phase variations may be disentangled to some extent. 



An evident property of the phase noise of J 17480 is that af- 
fects roughly three times more weakly the phases computed over 
the second harmonic than those of the fundamental. Together with 
the much easier description of the evolution of the second har- 
monic phases in terms of a physically plausible torque model (see 
Sec. 13.41 and in the following), this provides an observational in- 
dication of how the second harmonic component might be a bet- 
ter tracer of the spin frequency evolution of an accreting pul- 
sar. The enhanced stability of the second harmonic component 
with respect to the fundamental was already observed in a num- 
ber of AMSPs. Sudden jumps of the first harmonic phases not ac- 
companied by a corresponding variation of the second harmonic 
were observed during a couple of outburs ts of SAX J1808. 4-3658 
dBurderi et al.l l2006': Hartm an et ai1l2009l) . while during other out- 
bursts a comparable amoun t of no ise affected the two harmonic 
components dHartman et al.ll2008 a). A similar phenomenology is 
also observed from SWIFT JI749.4-2807 (Pa pitto et al., in prep.) . 
Further, in a coupl e of cases (XTE J18 07-295, lRiggio et al.l2008"bl : 
IGR J175 1 1-3057. 'Riggio et ai]|20I ll ) the noise affecting the sec- 
ond harmonic component was found to be significantly milder than 
that affecting the fundamental frequency; moreover, only by mod- 
elling the second harmonic component a physically plausible tim- 
ing solution could be proposed for these systems. While the dif- 
ferent stabi lity of the various harm onic components was taken into 
account by IHartman et al.l ( l2008ah to derive a minimum variance 
estimator for the phases of the pulse profile, not many attempts have 
been made to explain such a different behaviour of the phases com- 
puted over different Fourier components. If the observed variability 
of the pulse amplitudes and phases is due to irregular motion of the 
spot on the NS surface (L09) similar shifts should be observed in 
all the Fourier components. This has been obser ved and interpreted 
as such only in the case of XTE JI8I4-338 ( Papitto et all boOTi ; 
IWatts et alj|2008r) . While it could well be that substantial changes 
in the spot-shape or in the radiation-pattern of the emitted radia- 
tion determine a change of the pulse shape, this does not seem to 
naturally explain the enhanced stability of the second harmonic ob- 
served in a number of cases, and why only the use of the second 
harmonic phases as a tracker of the NS spin frequency generally 
produces results in accordance with theoretical expectations. 

A simple model proposed by iRiggio et al.l 1 120111) explains a 
similar behaviour in terms of modest variations of the relative in- 
tensity received by the two polar caps onto the NS surface, if the 
geometrical properties of the system make the signal coming from 
the two spots of similar amplitude. It is assumed that the pulse pro- 
file observed from an accreting pulsar is given by the sum of two 
signals with a similar harmonic content, emitted by two nearly an- 
tipodal spots on the NS surface (AA ~ vr, A7 ~ 27, where 7 and 
A are the latitude and longitude of the primary spot, respectively, 
and AA and A7 are the differences between the coordinates of the 
two spots), as it is the case if the magnetic field has a dipolar shape. 
The sum of the two signals (the total profile) will be the result of 
a destructive interference for what concerns the fundamental har- 
monic of the signal, since it is the sum of two signals with a phase 
difference of ~ tt. A constructive interference develops instead for 
the second harmonic of the total profile, since it is the sum of two 
signals with the same phase. The destructive interference regarding 
the fundamental frequency leads to (i) small values of the pulsed 
fraction of the first harmonic component, possibly comparable to 
that of the second harmonic, and (ii) large swings of the phase of 
the fundamental of the total profile due to modest variations of the 
relative intensity of the signals emitted by the two caps. The magni- 
tude of these effects is the largest when the spots on the NS surface 
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are exactly antipodal and produce signals that are observed at a 
similar amplitude. In this case swings up to 0.5 phase cycles can 
be shown by the phase computed over the fundamental frequency 
of the observed profile, without coiTespondingly large variations of 
the phase of the second harmonic. The maximum variability of the 
phase of the first haiTnonic component, as well as the maximum 
reduction of its amplitude, is obtained when the intensity of the 
signals observed from the two magnetic caps is very similar. This 
happens if the system is viewed nearly edge on and/or if the ro- 
tator is nearly orthogonal, that is with an angle between the spin 
and the magnetic axes of ~ 90°. Under these assumptions such a 
simple model can reproduce phase movements of the order of those 
observed from J 17480, as well as the enhanced stability observed 
from the second harmonic phases, in terms of variations of few per 
cent of the ratio between the intensity observed by the two spots. A 
thorough quantitative discussion of how this model applies to this, 
and other cases is deferred to a dedicated paper (Riggio et al., in 
prep.). 

In addition to the smaller residuals with respect to phe- 
nomenological models, an enhanced stability of the second har- 
monic phases is also indicated by the possibility of describing their 
behaviour by using a physically plausible model. In Sec. 13.31 and 
I3.4l we showed how the description of the pulse phases determined 
using the second harmonic component is significantly improved by 
assuming that the frequency derivative of the signal depends on a 
power close to one of the X-ray flux. The descritption of the pulse 
phases is further improved by th e assumption of a phase-flux cor- 
relation such as that proposed bv lPatruno et al .1 (2009). even if by a 
much lesser extent. On the other hand, the spin frequency derivative 
evaluated by considering the first harmonic phases shows oscilla- 
tions that are not correlated with the X-ray luminosity. 

Under the hypothesis that the second harmonic phases are a 
good tracker of the NS spin, we can infer some of the properties on 
the size of the inner radius of the accretion disc and on its depen- 
dence on the X-ray luminosity. To this end we consider the results 
presented in Sec. 13.41 where the phases of the second harmonic 
were fitted by assumung a relation, u ex L^ , for the dependency of 
the spin-up torque on the accretion luminosity. According to stan- 
dard accretion disc theory, the angular momemtum is transferred 
from the disc matter to the magnetosphere at the inner disc radius 
Rin. In turn, the size of the inner edge of the disc is determined 
by evaluating where the stresses exerted by the magnetic field lines 
on the accreting plasma become dominant with respect to the vis- 
cous stresses, that grant angular momentu m redistributi on in the 
disc far from the magnetosphere (see, e.g., lGhoshll2007L and ref- 
erences therein). The location of the inner disc radius is usually 
expressed as a fraction ^ of the Alfven radius evaluated in the case 
of spherical symmetric accretion 



„ „T f. 2/7 1,^1/7 0-2/7 r-2/7 4/7, 

Rin = 71 S.ri M, 1.4 -R,.io -^^37 ^27 km- 



(8) 



Here, ^27 is the magnetic-dipole moment in units of 10^^ G cm^, 
and the factor ^ is determined by capability of the field lines to 
thread the disc matter, and by the magnetic pitch, which in turn is 
limited b y the reconn ection-mechanisms of the field lines through 
the disc ( Ghosh & La mb 1979; Wang 1996, who quote values of ^ 
between 0.5 and I, respectively; se e also the discussion of these 
models given bv IBozzo et al.l 120091) . Inserting this expression in 
Eq. [7] it can be seen that the spin frequency derivative is expected 
to scale as !> oc Lj[ , if the only relevant torque is that exerted by 
matter orbiting in a Keplerian disc. 

The dependence on the X-ray luminosity of the spin frequency 



derivative evaluated from the phases of the second harmonic shown 
by J17480, (3 = 1.07 ± 0.03, is steeper than what is predicted 
by standard accret ion theory. Values of /3 close to unity were al- 
ready reported by Isildsten et al.l (119971) for the 0.47 s accreting 
pulsar in a LMXB, GRO J1744-28 (P = 0.957 ± 0.026) and 
the 105 s high field pulsar, A() 535-l-26 (J3 = 0.951 ± 0.026), as 
well as by IPaimar et al.l ( 119891) for the 42 s high field accreting 
pulsar, EXO 2030-1-375 (/3 = 1.21 ± 0.13). Expressing the de- 
pendence of the inner disc radius on the accretion luminosity as 
Rin oc L^", the estimate of j3 found for J 17480 translates in a 
value of Q = 2(1 — /3) = —0.14 ± 0.06. A negative value of 
Q would imply a positive correlation of the inner disc radius with 
the X-ray luminosity, which seems unlikely on physical grounds. 
On the other hand, a value close to zero indicates a very weak, if 
any, dependence of the inner disc radius on the accretion luminos- 
ity. A similar result, /3 = 0.97(4) [a = 0.06(8)], is obtained when 
adding to the model used to fit the phases a term depending linearly 
on the X-ray luminosity. 

By considering the value of the spin up torque, !>(f) = 
2.755(5) X 10^^^ Hz s^^, attained when the X-ray luminosity 
reached its peak value, Lx(f) = 9.23(4) x 10^'^ d^g erg s"\ 
and assuming that the NS accretes the Keplerian angular momen- 
tum of matter at the inner disc boundary Rin (Eq.[7](, the inner edge 
of the disc is estimated as _Ri„ (i) — 65(l)r;^jj_45 A/*,i.4i?76'^5^9 
km. A slightly smaller estimate, ~ 58 km, is obtained if the 
value of the frequency derivative obtained after adding a linear 
term in the X-ray luminosity to the model used to fit the phases, 
!>(f) = 2.602(6) X 10"^^ Hz s"\ is considered. The error on 
the distance (~ 8% of its estimate) is the largest source of un- 
certainty on the disc size measure. By taking it into account, the 
estimate of Rin{i) lies between 47 and 93 km (la confidence 
level). Inserting the measured values of Rin and Lmax in Eq. [8] 
yields to an estimate of the magnetic-dipole moment, fi27 = 

2.6(2) C''^^ d-^9 ll% K,w ^t[iA J?""^*. which has a simi- 
larly large uncertainty and for ^ ~ 0.5 varies between 5 and 16. A 
magnetic surface flux density between ~ 5 x 10^ and ~ 1.5 x 10^" 
G is therefore likely. 

Since JI7480 is a relatively slow pulsar, a coupling between 
magnetic field lines and plasma rotating in the accretion disc be- 
yond the inner disc radius transfers additional angular momentum 
to the NS. To take into account this effect, the torque on the NS can 
be expressed as 



N ^ 1.4 iVo 



- ujs/ujc 

1 — UJS 



(9) 



( lGhosh&Lamblll979l) . Here, uus = 2ttu^/rIJGM is the fast- 
ness parameter while uc is the critical fastness, which takes values 
between 0.35 and 1 depending on the description of t he interaction 
between the field and the disc ( JGhosh & Lamblll979l;IW ang 199^. 
Taking the value of ujc = 0.6 obtained bv lRomanova et alj l2002h 
from numerical simulations, and equating the observed torque to 
Eq.|9l the estimate of the inner disc radius at f reduces to roughly 
half of the estimate given above, Rin — 36 d^g km, where only 
the steep dependence on the distance was made explicit. The mod- 
elling of the phases of the second harmonic with a phisically plau- 
sible torque model yields therefore a val ue of the inner dis c ra- 
dius compatible with the estimate given bv lMiller et al.l ( 1201 ih on a 
spectroscopic basis, 41(5) Mt^iA km; the latter value was obtained 
by interpreting in terms of relativistic broadening the profile of the 
iron K-a line detected in the spectrum of the source observed by 
Chandra on 55493.445, (13.632 d since To), when the X-ray lu- 
minosity was a factor ~ 0.6 lower than at the peak. Still, these 
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estimates appear slightly larger than the inner disc radius implied 
by the observation of a QPO at 815 Hz during the outburst peak 
JAltamirano et al.ll201(]n . Rin — 19 M«,i.4 km, if such a feature 
is interpreted in terms of a Keplerian orbital frequency. Given the 
steep dependence of our inner disc radius estimate on the distance 
to the source, the size given by the Keplerian interpretation of the 
QPO would be reconciled with the value indicated by the torque 
analysis presented here if a distance of ;i, 7 kpc to the Globular 
cluster is assumed. A similar value is marginally compatible with 
the most recent determination (5.9 ± 0.5 kpc), though it has to be 
noted that heavy, differential reddening affecting observations of 
this Globular clus ter make the determi nation of its distance some- 
what problematic JLanzoni et al.l2010l and references therein). 
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Table 3: Best-fitting parameters of the spectral continuum evaluated from the 2.5-30 keV spectra observed by PCA on-board 
RXTE. Columns are: Observation ID (labels refer to the groups of observations defined in text); Start time of the observation 
(MJD since To = MJD 55479.81229); End time of the observation (MJD since To); Exposure (ks); kTe (keV); t; kTs (keV); 
Unabsorbed flux in the 2.5-25 keV band (10~^ erg cm~^ s~^); Bolometric luminosity calculated by extrapolating the best- 
fitting model to the 0.05-150 keV energy band (di.g 10^^ erg cm~^ s~^); Reduced chi-squared of the fit. The quoted errors are 
evaluated at the 90% confidence level on each parameter (Ax^ = 2.706). 
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Table 3 - continued from previous page 
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Table 4: Best-fitting parameters of the spectral continuum evaluated from the combined spectra composed of 
(quasi)simultaneous observations performed by XRT on-board Swift (1-10 keV) and by PCA on-board RXTE (2.5-30 keV). 
Columns are: Observation ID; Start time of the observation (MJD since To = MJD 55479.81229); End time of the observation 
(MJD since To); Exposure (ks); A''^ (10^^ cm"^); fcTe (keV); r; kTs (keV); XRT/PCA normalization; Unabsorbed flux in the 
2.5-25 keV (PCA) and in the 1-10 keV (XRT) band (10~* erg cm~^ s~^); Bolometric luminosity calculated by extrapolating 
the best-fitting model to the 0.05-150 keV energy band (dl g 10^'^ erg cm~^ s^^). The quoted errors are evaluated at the 90% 
confidence level on each parameter (Ax^ = 2.706). 
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Table 5: Best-fitting parameters of the spectral continuum evaluated by combining the 1-10 keV spectrum observed by XRT 
on-board Swift starting on MJD 55479.801, with the 4.5-25 keV JEM-X spectrum and the 20-150 keV ISGRI spectrum ac- 
cumulated during the INTEGRAL revolution No. 975 (MJD 55479.365 - 55479.548). This group of observations is labeled 
in text as Obs. [A]. Columns are: Observation ID or satellite revolution; Start time of the observation (MJD since Tq = MJD 
55479.81229); End time of the observation (MJD since To); Exposure (ks); N^ (10^^ cm"^); kT^ (keV); r; kTs (keV); 
Unabsorbed flux in the 1-10 keV (XRT), 4.5-20 keV (JEM-X) or in the 20-150 keV (ISGRI) bands (10~* erg cm"^ s"^); 
Normalisation of JEM-X and of the ISGRI spectra with respect to XRT; Reduced chi squared of the fit. The quoted errors are 
evaluated at the 90% confidence level on each parameter (Ax^ ~ 2.706). 



ObsID (Sat. rev.) Tstart (MJD) T^nd (MJD) Expo (ks) iV^ (10^^ cm~^) kTe (keV) 



kTg (keV) 



(JEM-X - ISGRI)/XRT x /dof 



00031841002 


-0.0103 


0.0127 


2.0 


975(JEMX) 


-0.4471 


-0.2640 


8.7 


975(ISGRI) 


-0.4477 


-0.2821 


8.2 



1.1(2) 



2o: 



4.4 



TTTT- 



0.7(1) 



0.053(6) 
0.05(1) 
0.08(1) 



109/109 



0.68 



+0.17 
0.15 



Table 6: Best-fitting parameters of the spectral continuum evaluated from the combined 4.5-150 keV spectra observed by JEM- 
X and ISGRI on-board INTEGRAL. Columns are: Satellite revolution; Start time of the observation (MJD since To = MJD 
55479.81229); End time of the observation (MJD since To); Exposure (ks); kTe (keV); r; kTs (keV); Unabsorbed flux in the 
4.5-20 keV (JEM-X) or in the 20-150 keV (ISGRI) band (10^® erg cm"^ s"^); ISGRI "soft" unabsorbed flux in the 20-50 
keV band (10"** erg cm"^ s"^); ISGRI "hard' unabsorbed flux in the 50-150 keV band (10"* erg cm"^ s"^); Reduced chi 
squared of the fit. The quoted errors are evaluated at the 90% confidence level on each parameter (Ax^ ~ 2.706). 
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The quoted fluxes are evaluated by fitting the ISGRI points with a power law model. 
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